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Abstract We investigate Coulomb blockade thermometers (CBT) in an intermediate tem- 
perature regime, where measurements with enhanced accuracy are possible due to the in- 
creased magnitude of the differential conductance dip. Previous theoretical results show 
that corrections to the half width and to the depth of the measured conductance dip of a 
sensor are needed, when leaving the regime of weak Coulomb blockade towards lower tem- 
peratures. In the present work, we demonstrate experimentally that the temperature range of 
a CBT sensor can be extended by employing these corrections without compromising the 
primary nature or the accuracy of the thermometer. 
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1 Introduction 

Since the single electron tunneling phenomena were described about 30 years ago 1 1^, many 
devices based on small tunnel junctions have been demonstrated. An example of such a 
device is the Coulomb blockade thermometer (CBT) using one |2| or two dimensional ar- 
rays of tunnel junctions |3jj4j. CBT based on electron tunneling enables one to relate the 
measured voltage to temperature with the help of Boltzmann constant kg, which could thus 
provide a means to revise the definition of kelvin ||5]- A related idea forms the basis of 
shot-noise thermometry, described in Refs. ||6]|7]. 
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2 Theoretical background 

CBT thermometry is based on the change of electrical conductance G of tunnel junction 
arrays. Typically, CBT works in a weak Coulomb blockade regime Ec kgT , where the 
charging energy of the system with A' junctions in series is Ec = [{N — 1) /N]e^ /Cz and Cz 
is the total capacitance of an island between the junctions that depends on the physical size 
of the contact and the self capacitance of the island. In this regime, CBT sensor works as a 
primary thermometer, meaning that calibration is not needed |2|. For a uniform array, the 
differential conductance scaled by its asymptotic value Gt at large positive and negative 
voltages can be written in terms of bias voltage V as 



where the function g{x) = [jcsinh(jt:) — 4sinh^(x/2)]/[8sinh^(x/2)] determines the bias de- 
pendence. 

Primary measurements of temperature can be achieved from the measurement of full 
width of the conductance dip at half minimum 



Thus, obtained temperature does not depend on the geometry or material of a sensor. A 
secondary temperature measurement can be obtained by recording the depth of the zero bias 
conductance dip 



with parameter mai = Ec/ksT . As was shown in Ref. [8], at high temperature (in that case 
4.2 K) the measured conductances follow quite well the expression in Eq. ((Sjl. 

In this letter, we demonstrate experimentally the extension of Coulomb blockade ther- 
mometry into an intermediate temperature regime described theoretically in Refs. Il8]|9l- Ex- 
periments employing quantum dot thermometry II101I11II12I[T31 . investigating electrical con- 
duction in metallic nanodots 1141 or granular films III5I work in this intermediate Coulomb 
blockade regime as well. 

The intermediate regime is defined by Ec ~ fc^T, between the traditional CBT operation 
range Ec ^ keT, and the full Coulomb blockade Ec ^ kgT, see Fig. [T] Solid black and 
dash - dotted blue lines correspond to the numerically calculated temperature dependent 
conductance of CBT at zero bias using theory of single electron tunneling |2| (refered to 
"full theory" in the following text). The two limiting cases, "minimum blockade" (dash - 
dotted blue line) and "maximum blockade" (solid black line), represent the two opposite 
extremes of the possible background charge configurations. 

Thus we compare the two extreme cases to indicate the uncertainty due to the unknown 
background charge values. As long as the two curves follow each other, the uncertainty due 
to background charges can be neglected. In the regime Ec ^ ksT the charge sensitivity plays 
a dominant role |i61, as the measured conductance depends on the randomly fluctuating 
background charge configuration, see Fig. [T] To fully tune the background charges in a large 
array one should have an equally large number of gate controls. Instead of using elaborate 
gate control, we therefore limit the use of the thermometer to a regime, where the influence 
of background charges on temperature reading stays sufficiently small. 

The intermediate temperature regime is beneficial for thermometry due to enhanced 
signal magnitude <x Ec/ksT . In addition, the harmful contribution of nonlinear I-V curves 




(1) 



Vii2'^5A'i9NkBT /e. 



(2) 



Primary thermometry in the intermediate Coulomb blockade regime 



3 





■ 12 


o 


9 


^0.1 




< 






■ 




0.01 - 



SET 



maximum blockade 

— minimum blockade 

— - "j^ 

— • uJ6-u^''/60+u^'/630 



\ 



CBT 



L 



iJ- 



J I 



0.03 0.1 



k T/E 

B c 



10 



Fig. 1 (Color on-line) Theoretical temperature dependence of the depth of the zero bias conductance dip of 
a tunnel junction array influenced by Coulomb interaction. The left side of the plot represents the working 
regime of single electron devices (SET), where E,- S> ksT. Colored area represents an intermediate regime, 
where hgT ^ E^, which is the focus of this work. On the right side of the plot the normal working regime 
of CBT sensors is represented, where E^ <S ksT. Two extremes of the possible experimental results due to 
unknown background charge configuration are represented by maximum (solid black line) and minimum 
(dash - dotted blue line) blockade curves. The third order result, Eq. ^5), is presented by the dash double - 
dotted green line. The inset represents the error in temperature reading which one makes using the third order 
approximation of Eq. (5). The orange dotted curve represents the relative difference between the numerically 
calculated minimum blockade (dash - dotted blue line, Fig.[T) curve and the third order curve (dash double - 
dotted green line, Fig.^. The sohd green line on the inset plot represents the relative difference between the 
maximum blockade (solid black line) and the third order approximation curve. 



(predominantly due to finite tunnel barrier height) as a function of applied bias voltages is 
less abundant |17|. Another motivation of the present work is to expand the temperature 
range where a single CBT sensor can be used [ 18 1. 

To obtain the temperature with high accuracy, one has to make a series expansion up to 
higher orders in towards low temperatures [8 1. In this work, we show experimentally that 
the temperature of a CBT sensor can be obtained by employing corrections to Eqs. (|2} and 
^ as 

Vu2 = 5.439NkBT{l +03921^) /e, (4) 

(jT 

1 1 2 1 3 

Note that as ¥1^2 ™d AG/Gj are measured quantities, Eq. (|4} maintains the primary na- 
ture Q. We examine here the series expansion up to the third order, which gives an accuracy 
of 1.1 % (see in the following text). Since we assume that the scatter in the resistance of the 
turmel junctions can bring a comparable eiTor, we do not consider higher order contributions. 



' The prefactor (0.3921) in front of the depth of zero bias conductance dip, AG/Gt in Eq. (4), is the 
original value and the same as in Ref. 181. In Ref. 191 there is a typo (0.3992). 
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The solution of Eq. ([5} in terms of conductance is given by 
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In Fig.[TJ one can see that Eq. ^ (dashed red Une) at low kgT /Ec is no longer valid. 
The exact dependence of AG/Gj is shown by solid black and dash - dotted blue lines. On 
the other hand, the higher order series expansion presented by Eq. ([5]l (dash double - dotted 
green line) follows the exact dependence with good accuracy down to lower temperatures. 
In the inset of Fig. [T] the orange dotted curve represents the relative difference between 
minimum blockade and the third order approximation expressed in Eq. ([S}. The solid green 
curve represents the relative difference between maximum blockade and the third order ap- 
proximation. To keep the uncertainty in the measurement small, we limit the range of CBT 
operation to above ksT /Ec ~ 0.4. This means that the third order result deviates less thein 
2.5% from either of the opposite extremes (maximum and minimum blockade) of possible 
experimental results. Note that although the expansion parameter m^? assumes values ex- 
ceeding unity, the series of Eq. (|5]l does not diverge even in this case because the multipliers 
diminish quickly for higher orders 0. 



3 Experimental Details 

The measured sensors consist of parallel arrays of tunnel junctions in series. Figure |2] (a) 
shows an image of three junctions. 




Fig. 2 (Color on-line) The tunnel junctions and the measurement setup, (a) In an optical microscope image 
of one of the sensors, three tunnel junctions (circles) formed by overlapping area of two, top and bottom, Al 
fingers (wide orange lines) are shown, (b) Close-up SEM image of a cross-section of a tunnel junction. The 
thin AI2O3 layer between the two 250 nm-thick Al electrodes is the tunnel banier The tunnel junction has 
round shape with a diameter of 500 nm. (c) Schematic diagram of the measurement setup. A CBT sensor 
consists of 9 parallel an'ays of 49 tunnel junctions in series, which connect to the bonding pads. On the left 
side of a schematic diagram is the sensor, which is placed in the refrigerator The right side of a schematic 
diagram represents the measurement setup at room temperature. 



^ A more natural expansion parameter would be the depth in the lowest order Wfj = mjv/6. In this case 
Eq. {5j would read AG/Gt = wn- (3/5)w^ + (12/35)w^ 
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These sensors were made by ex-situ tunnel junction fabrication method, which utilizes 
plasma etched via in a dielectric (SiOx) for defining the tunnel junctions. Fabrication process 
is described in detail in Ref. 1 19|. Tunnel junctions were configured in a round shape of 500 
nm in diameter (see Fig.|2](b)). In order to keep the size of the tunnel junctions identical for 
sensors with different Ec, the charging energy of a sensor was tuned by varying the overlap 
area of the cooling fingers. In Fig.|2](a), for the sensor with the highest charging energy the 
overlap area is very small and noticeable only around the junction area. For the sensor with 
the lowest charging energy, the top and the bottom cooling fingers are totally overlapping 
(not shown). 

The sensors were measured in a ^He — ^ He dilution refrigerator. Detailed schematic 
diagram of the measurement setup is presented in Fig. |2](c). Differential conductance was 
measured, as shown in Fig. |2] (c), using SR 830 Lock-In amplifiers and DL instruments 
low noise voltage (1201) and current (1211) preamplifiers. In addition to voltage biasing 
V = VgenRi/iRi +Ri), a small (2% of V1/2) voltage was applied. To obtain the differential 
conductance G = 51/ 5V, both the ac voltage 5V and current 51 were measured. 

4 Results and Discussion 

In the measurement we focus on the intermediate regime, keT ~ Ec- We performed the 
measurements for three CBT sensors with different values of Ec = 45 mK (sensor A), 
Ec = 94 mK (sensor B) and Ec = 129 mK (sensor C) in a temperature range between 
55 mK and 470 mK . In case of sensor C we used a low-charging energy CBT sensor as 
a reference temperature thermometer (sensor REF, Ec = 22 mK). It has a tunnel junction 
size of 800 nm in diameter and was operating in the linear CBT working regime, where the 
second and third order eifects are small. 

In Fig. [3] (a), we present the measured normalized conductance of sensors A and B vs 
applied voltage at a refrigerator temperature of 55 mK. The experiment shows the expected 
difference of the conductance peak depth corresponding to the two charging energies. 

In Fig. [3] (b), we present the experimental conductance curve (blue dots) for sensor B 
and the full theory fit to the data (dashed red line). The theoretical model was calculated as 
follows: first, I — V curves were obtained by employing the master equation, as shown in 
Ref. O. In addition, we take overheating effects into account ||9l- The electron temperature 
is higher than the lattice temperature, due to the weakness of the electron-phonon interaction. 
These effects are relatively small due to the cooling fingers (see Fig. [2] (a)), which provide 
a good thermalization between the arrays and the substrate [j. Finally, we differentiate the 
numerical I — V curves to fit them to the measured conductance curves. 

The charging energies and the temperatures were obtained using a standard technique: 
the sensors were measured at higher refrigerator temperatures (between 420 - 470 mK) and 
fitted to numerically calculated full conductance curves, as described above. At low tempera- 
tures, assuming constant charging energies, we extract the temperatures of the sensors using 
higher order series expansion expressed in Eq. (|4]. During our measurements, we record the 
temperatures of the fridge via a calibrated Ru02 thermometer. At any temperature we have 
the possibility to compare the Ru02 value and the temperature given by the fitting procedure 
for our sensor. 

The values of Ec = 94 mK and T =47.1 mK obtained from the full theory fit to sensor B 
data in Fig.[3](b) were used to calculate the depth of the conductance dip for the first (Eq. (|3}) 

^ At ±^1^2/2, the bias points used in the determination of the temperature, the islands are overheated by 
about 8 mK at maximum. 
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Fig. 3 (Color on-line) Measured conductance of the CBT sensors at 7" = 55 mK detemiined by R11O2 and 
the theoretical predictions, (a) Normalized conductances G/Gt vs applied voltage V measured for sensors A 
and B. (b) The experimental data for sensor B (blue dots) together with the numerically calculated full theory 
fit on top (dashed red line). Minimum of the measured conductance at zero bias, which also corresponds to 
the minimum of the numerically calculated full theory curve, marked as red square, and the minimum of the 
conductance dip calculated using the third order approximation marked as blue star. The solid green line of 
Eq. Q corresponds to the first order curve and the rhombus to its minimum. The inset shows the full theory 
curve (red) fitted to measured conductance curve (blue dots) and the first order curve (solid green line) in the 
bias range near ¥1/2/'^- 

and the third (Eq. (|5]l) order approximations. The experimental data yield ZiG/Gj- = 0.2815, 
whereas the third order approximation utilizing Eq. (|5) gives AG/Gt = 0.2789. These val- 
ues correspond to the red square and the blue star in Fig. |3](b), respectively. As one can see 
in Fig.[3](b), the numerically calculated full theory curve follows the experimental conduc- 
tance curve very well, and thus our experiment is in agreement with the theoretical model 
of single electron tunneling. The value of the depth of the experimental conductance curve 
(red square) also corresponds to the depth calculated using the third order approximation 
(blue star), Eq. (|5]l. The series approximation differs from the result of the full theory by 
1.1% only. The value for the depth of the conductance dip, using only the first order series 
expansion expressed in Eq. is AG/Gt = 0.3326. To show the difference between the 
first and the third order series expansions, we plot the first order curve (Eq. ([T)) on the same 
graph. The shape of the first order conductance curve is significantly different from the ex- 
perimental (blue dots) and theoretical (dashed red line) curves. This difference would lead 
to an error of about 30% in temperature reading, if the first order result was employed. 

In Fig. |4] (a), we present the first, second and third order series expansions (dashed 
black, dash dotted blue and solid red curves, respectively), calculated using Eq. (|5]l. On the 
same plot we show our experimental data points for the depth of the conductance dip for 




V/N [^V] 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 11 






















#0.89 










• 


experiment K 


M 0.88 












full theory \ 


/ 0.87 












1st order 1 










■ 


AG/G^ experiment % 


f 0.86 










♦ 


AG/G.^ 1st order W 


f 0.85 










★ 


AG/G.^3rd order M 
1 1 1 1 1 1 1 1 1 T 


^ 0.84 
, , . 1 








1 




22 24 
1 1 1 


26 28 30 
1 , . . 1 



Primary thermometry in the intermediate Coulomb blockade regime 



7 



sensor B. The temperature axis is given by the full theory fit to the measured conductance 
curve. As one can see, the experimental data fit to the third order curve very well at low 
temperatures, when Ec = ksT . In Fig.|4](b), we present the depth of the conductance curve 
vs. temperature for sensors A, B and C. We show the experimental data, marked as circles, 
stars and rhombi, together with the first (dashed black line) and the third (solid red line) 
order series expansions in 40 - 80 mK temperature range. The data follow the third order 
curve, calculated using Eq. ([5]l. 




T [mK] 

Fig. 4 (a) The depth of the conductance dip vs temperature for sample B. The three Unes correspond to the 
first (dashed line), the second (dash - doted line) and the third (solid line) order series expansions, (b) Close- 
up of (a), zoomed to the low temperature regime, for three CBT sensors. From top to bottom the sets of curves 
conespond to sensors C, B and A, respectively. The dashed black line corresponds to the first order and solid 
red line to the third order series expansion. Experimental data are shown as circles, stars and rhombi. 

In Fig. |4] (b), the top two curves together with experimental data correspond to sensor 
C. We provide an additional proof realizing the temperature scale independently using low- 
charging energy CBT sensor as reference. Sensor REF, operating in the traditional CBT 
range, Ec <S kgT , was measured at the same time with sensor C. For the latter sensor data 
points also follow the third order curve. The temperatures of the two sensors agree within 
0.5%. Sensor REF had been earlier compared favorably with the PLTS 2000 temperature 
scale f20| with an accuracy of 1%. The difference in depths of the conductance dips between 
the first and the third order series expansions for sensor REF is less than 4%. 

The middle and the bottom sets of curves in Fig.|4](b) correspond to the sensors B and 
A, respectively. The temperature axis for the two bottom sets of experimental data of these 
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sensors (stars and rhombi) was determined by the full theory fit to our data (as before in 
Fig. 11(a)). In this case, the fitting forces the depth of the conductance dip to be in accor- 
dance to the temperature. We do not observe any deviations of the fitted curve within the 
experimental accuracy (see inset of Fig. |3](b)). This nearly perfect agreement leaves only 
a small deviation between the measured depth of the conductance dip and Eq. ([5]l. The de- 
viation lies within the uncertainty of the value of Ec, especially when its absolute value 
becomes smaller, as with sensor A. 



5 Summary 

To summarize, we demonstrate primary Coulomb blockade thermometry operated in a wide 
range extending to the Ec ~ kgT regime, where the observable signal magnitude and con- 
sequently the accuracy are enhanced. We conclude that it is possible to operate the CBT 
thermometer reliably down to kgT /Ec = 0.4, by employing the third order series expan- 
sion. Finally, the accuracy of this method is limited by the influence of background charge 
distribution, but stays below 2.5% down to kgT /Ec = 0.4. 
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